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PHENIX has measured the centrahty dependence of mid-rapidity pion, kaon and proton transverse 
momentum distributions in d+Au and p+p coUisions at ^/snn = 200 GeV. The p+p data provide a 
reference for nuclear effects in d+Au and previously measured Au+Au collisions. Hadron production 
is enhanced in d+Au, relative to independent nuclcon-nucleon scattering, as was observed in lower 
energy collisions. The nuclear modification factor for (anti) protons is larger than that for pions. 
The difference increases with centrality, but is not sufficient to account for the abundance of baryon 
production observed in central Au+Au collisions at RHIC. The centrality dependence in d+Au 
shows that the nuclear modification factor increases gradually with the number of collisions suffered 
by each participant nucleon. We also present comparisons with lower energy data as well as with 
parton recombination and other theoretical models of nuclear effects on particle production. 
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I. INTRODUCTION 

Since the early 1970's, it is well established P, 0, 
0], that energetic particle production in proton-nucleus 
(p+A) collisions increases faster than the number of bi- 
nary nucleon-nucleon collisions. This effect, called the 
"Cronin effect" , is a manifestation of the fact that parti- 
cle production and propagation is influenced by the nu- 
cleus. If the A-dependence of the invariant cross section, 
/, of particle i in p+A collisions is parameterized as 

hipT,A)=h{pT,l)-A''^^PT) (1) 

then it has been observed that at is greater than 
unity above some transverse momentum value, typi- 
cally 1-1.5 GeV/c, denoting significant enhancement of 
particle production in p+A collisions. The enhance- 
ment depends on the momentum and the type of par- 
ticle produced, with protons and antiprotons exhibit- 
ing a much larger enhancement than pions and kaons 
at > 2-3 GeV/c. At y^s^N = 27.4 GeV, 
the enhancement peaks at around pt=4.5 GeV/c, with 
aji+ ~ ajr+ = 1.109 ± 0.007, while, at the same momen- 
tum, the protons can be described by an a-factor of 
ap-a^+ =0.231 + 0.013 0. 

Although the observables in Eq. ^ have been clearly 
related to the nuclear medium, the cause of the Cronin 
enhancement and its species dependence are not yet com- 
pletely understood and further experimental study is 
warranted in its own right. Furthermore, in the search 
for the Quark Gluon Plasma at the Relativistic Heavy 
Ion Collider (RHIC), the Cronin effect is extremely im- 
portant, as novel effects observed in central Au+Au colli- 
sions require good control of the initial state conditions. 
At RHIC energies, it was discovered that hadron pro- 
duction at high transverse momentum {pT > 2GeV/c) 
is suppressed in central Au+Au collisions ^ compared 
to nucleon-nucleon collisions. Such suppression may be 
interpreted as a consequence of the energy loss suffered 
by the hard-scattered partons as they propagate through 
the hot and dense medium. However, since the Cronin 
effect acts in the opposite direction enhancing the hadron 
yields, it has to be taken into account when the parton 
energy loss is determined from the data. 

Another discovery at RHIC, originally unexpected, is 
that the yields of p and p at intermediatepr (1.5 < pt < 
5GeV/c ) in central Au+Au coUisions are com- 
parable to the yield of pions, in striking contrast to the 
proton to pion (p/tt) ratios of ~ 0.1 - 0.3 measured in p+p 
collisions [g]- Novel mechanisms of particle production 
in the environment of dense matter, such as recombina- 
tion of boosted quarks Q or contributions from baryon 
junctions 0, which can become dominant in the pres- 
ence of pion suppression were proposed to explain the 
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data. Since it has been observed that at lower energies 
Cronin enhancement is stronger for protons than for pi- 
ons 0, this effect has to be considered at RHIC before 
new physics is invoked. 

The effects from the initial state are best studied by 
performing a control experiment in which no hot and 
dense matter is produced. Deuteron + gold collisions 
at y/SNN~ 200 GeV serve this purpose. Since there is 
no hot and dense final state medium, the initial state 
conditions become accessible to the experiment. In addi- 
tion to Cronin enhancement, known initial state effects 
also include nuclear shadowing and gluon saturation [ij . 
The Cronin enhancement is usually attributed to mo- 
mentum broadening due to multiple initial state soft [l^ 
or semi-hard 0, Q, 0, 0| scattering. Such models 
typically do not predict the particle species dependence 
observed in the data. Recently, Hwa and collaborators 
provided an alternative explanation due to final state in- 
teractions. The particle species dependent enhancement 
is attributed to recombination of shower quarks with 
those from the medium, where no distinction is made 
if hot or cold nuclear matter is produced ,17j . Identified 
hadron production measured as a function of centrality 
brings important experimental data relevant to the long 
outstanding problem of the baryon Cronin effect. The 
dependence of the enhancement upon the thickness of 
the medium, or the number of collisions suffered by each 
participating nucleon, can help differentiate among the 
different scattering models, and the species dependence 
helps to separate initial from final state effects in d+Au. 

The paper is arranged as follows. Section II describes 
the experiment, data analysis, and systematic uncertain- 
ties. Section HI presents hadron spectra, yields and the 
resulting nuclear modification factors. Discussion of the 
centrality, energy and species dependence of the nuclear 
modification factors and implications for understanding 
of the Cronin effect are in section IV. Section V presents 
conclusions. 



II. EXPERIMENT AND DATA ANALYSIS 
A. Data Sets and Trigger 

Data presented here include collisions at y/sJjN ~ 
200 GeV of Au+Au taken in the 2002 run of RHIC and 
d+Au and p+p collected in 2003. In the following we 
discuss analysis of the p+p and d+Au data; details of 
the Au+Au analysis, and the Au+Au results, are found 
in 7J. Events with vertex position along the beam axis 
within \z\ <30 cm were triggered by the Beam-Beam 
Counters (BBC) located at |r;| = 3.0-3.9 T^. The min- 
imum bias trigger accepts 88.5 ± 4% of all d+Au colli- 
sions that satisfy the vertex condition, and 51.6 ± 9.8% 
of p+p collisions. A total of 42 x 10^ minimum bias 
d+Au events and 25 x 10^ minimum bias p+p events 
were analyzed. 

In p+p collisions, PHENIX determines the differential 
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TABLE I: Mean number of binary collisions, participating nucleons from the Au nucleus, number of collisions per participating 
deuteron nucleon, and trigger bias corrections for the d+Au ccntrality bins. 



00-20% 



20-40% 



40-60% 



60-88% 



(Npart) 
{V = Ncou/N^art) 

Trigger bias correction 



15.4±1.0 
15.6±0.9 
7.5±0.5 
0.95±0.3 



10.6±0.7 
11.1±0.6 
5.6±0.4 
0.99±0.007 



7.0±0.6 
7.7±0.4 
4.0±0.3 
1.03±0.009 



3.1±0.3 
4.2±0.3 
2.2±0.2 
1.04±0.027 
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C^BBC \_ 



■ciTj{PT)-cgg^ 
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The BBC cross section ub rc was determined via the 
van der Meer scan technique 19]. In this p+p data set, 
asBC = 23.0 ± 2.2(9.6%) nib. A^|^*g' is the total number 
of BBC triggers analyzed. The factor C^Jj°j{pT) denotes 
the efficiency and geometrical acceptance correction, cal- 
culated with a detailed GEANT Monte Carlo simulation 
113 of the PHENIX detector. Cf^^(pT) normalizes the 
cross section in one unit of rapidity and full azimuthal 
coverage. The C^^*^ factor corrects for the fact that the 
forward BBC trigger counters measure only a fraction 
of the inelastic p-|-p cross section. This subset of events 
on which the BBC triggers contains only a fraction of 
the inclusive particle yield at mid-rapidity. For charged 
hadrons this fraction was determined using triggers on 
the beam crossing clock; the fraction was found to be 
0.80 ± 0.02, independent of pr 
our nomenclature, the inverse of this fraction 

In d+Au collisions PHENIX measures the inelastic 
yield per BBC triggered event. ^ The collision ccntral- 
ity is selected in d-|-Au using the south (Au-going side) 
BBC (BBCS). We assume that the BBCS signal is pro- 
portional to the number of participating nucleons (N^^f) 
in the Au nucleus, and that the hits in the BBCS are un- 
correlated to each other. We use a Glauber model 
and simulation of the BBC to define 4 centrality classes 
in d+Au collisions, as discussed in detail in [23. The 
deuteron nucleus is modeled after a wave function due to 
Hulthen M 



The C^f^<^ term is, in 



aP{a + (3) 
27r(a - /3)2 



(e- 



)/^pn (3) 



where a = 0.228 fm~^ and (3 = 1.18 fm~^ and refers 
to the separation between the proton and the neutron. 
The Au nucleus is modeled using a Woods-Saxon density 



^ To convert the reported inelastic yield to differential cross sec- 
tion, one must multiply by the beam-be am trigger cross section 
"""(dAu) = 1.99 ± 0.10b 0. 



p{r) = 



Po 



1 



exp I 



(4) 



where a = 0.54 fm and c = 1.12 x A^/^ - 0.86 x A'^^^ = 
6.40 fm, a value which agrees well with the measured 
charge radius of i? = 6.38fm for gold. 

Using the above parameters and taking into account 
the BBC efficiency the mean number of binary collisions 
along with the mean number of participating nucleons 
from the Au nucleus that correspond to each centrality 
bin are shown in Table I. For the minimum bias d-f Au 



collisions, (Ncou) = 8.5 ± 0.4 and {Nj 



part/ 



9.1 ± 0.4. 



of a\%(dAu)e 



BBCf 
MB 



As for p-l-p collisions, there is a BBC trigger bias, but 
it is much smaller in d+Au. In addition, a second bias 
occurs in d+Au centrality selected collisions. This sec- 
ond bias arises from the fact that events containing high 
Pt hadrons from hard scatterings may have larger mul- 
tiplicity, and consequently produce a larger signal in the 
BBCS. Such events would be considered more central 
than events without a hard scattering. This effect gives 
an opposite bias from the first trigger bias effect in the 
most peripheral bin as events can be shifted out of this 
bin but not into it. We correct for both biases in d+Au 
collisions using simulations and a Glauber model. The 
combined corrections for these effects range from 0-5%, 
depending on the centrality category, and are shown in 
Table I. Systematic uncertainties on these corrections are 
less than 4%. 



B. Tracking and Particle Identification 

Charged particles are reconstructed using a drift cham- 
ber (DC) and two layers of multi-wire proportional cham- 
bers with pad readout (PCI, PCS) Pattern recogni- 
tion is based on a combinatorial Hough transform in the 
track bend plane, while the polar angle is determined by 
PCI and the location of the collision vertex along the 
beam direction 25]. The track reconstruction efficiency 
is approximately 98%. 

Particle momenta are measured with a resolution 
Sp/p = 0.7% e l.l%p (GeV/c). The momentum scale 
is known to 0.7%. Particle identification is based on par- 
ticle mass calculated from the measured momentum and 
the velocity obtained from the time-of-flight and path 
length along the trajectory. The measurement uses the 



5 




Time of Flight [ns] 



FIG. 1: (color online) Time of flight vs. charge/momentum 
for particles in minimum bias d+Au collisions. Bands corre- 
sponding to the different hadrons are labeled. Electron and 
muon bands are also visible, though not labeled. 
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FIG. 2: (color online) Mass squared distribution for posi- 
tively charged tracks with 1.2 < pr < 1.6 GeV/c, using the 
high resolution time-of- flight measurement. The solid regions 
indicate the mass squared ranges for accepted pious, kaons 
and protons, respectively, from left to right. 



portion of the spectrometer containing the high resolu- 
tion time-of- flight (TOF) detector, which covers pseudo- 
rapidity -0.35 < rj < 0.35 and Acjj = tt/8 in azimuthal 
angle. The timing uses the BBC for the global start, and 
stop signals from the TOF scintillators located at a radial 
distance of 5.06 m. Tracks are required to match a hit 
on the TOF within ± 3 standard deviations {a) of the 
track projection to the TOF radial location. The particle 
yields are corrected for losses due to this matching cut. 

Figure^illustrates the performance of the hadron iden- 
tification system. The system resolution is ti ~ 130 ps for 
both d+Au and p+p; the start time resolution is limited 
by the low hit multiplicity in the BBC. The TOF reso- 
lution achieved allows for a clear separation of tt/K and 
K/p up to pt — 2 GeV/c and px — 3.6 GeV/c, respec- 
tively. A 2(7 cut in mass squared is used to separate the 
different hadron species, as shown in Fig. |21for hadrons 
in the pT range from 1.2 - 1.6 GeV/c. A clear separation 
between pions and kaons is seen. The particle yields are 
corrected for losses due to the 2a cut. 

Corrections to the charged particle spectrum for geo- 
metrical acceptance, decays in flight, reconstruction effi- 
ciency, energy loss in detector material, and momentum 
resolution are determined using a single-particle GEANT 
Monte Carlo simulation. 

The proton and antiproton spectra are corrected for 
feed-down from weak decays via a Monte Carlo simula- 
tion using as input experimental data on A production. 
The total number of protons produced in the collisions 
can be written as: p+OM{A+T.°+E°+E- +n-)+0.52J:+ 
where p denotes the primordial number of protons pro- 
duced, and the other symbols denote the primordial num- 
ber of those particles produced in the collision. 0.64 
and 0.52 are the branching ratios for A pn" and 
E"*" — > p7r", respectively. The hyperons listed together 
with A decay to A with approximately 100% branching 



ratio, and so yield protons with the A —f pm branching ra- 
tio. We estimate the proton and antiproton spectra from 
weak decays and subtract them from the measured yields, 
using experimental data from UA5 ,.26j from non-single 
diffractive ^/s — 200 GeW p+p collisions, and preliminary 
A and A spectra from STAR in p+p and d-|-Au collisions 
at the same energy 27, 28, 29, 30] . The shape of the 
S", and S~ spectra are constructed from the A spec- 
trum by mj' scaling (i.e. under the assumption that these 
hadrons are all produced with roughly the same spectrum 
in transverse mass [soj]'). The relative normalization of 
k,YP,'EP, and S" from UA5 ^ is used for both p+p 
and d+Au collisions. This is justified by the similarity of 
the Cronin effect for different baryons. The contribution 
of f2~ is negligible, and is not included in the correction. 
The Monte Carlo simulation decays these baryons and 
propagates the products through the PHENIX magnetic 
field and central arm detectors, accounting also for the 
change of momentum distributions between parent and 
daughter particles due to decay kinematics. The result- 
ing proton and antiproton spectra are then subtracted 
from the measured inclusive spectra. The fractional con- 
tribution from feed-down protons from weak decays to 
the total measured proton spectrum, /, is approximately 
30% at moderate and high px growing slowly for lower 
Pt and reaching 40% at pT = 0.6 GeV/c. / is shown as 
a function of pT in Fig. (SJ along with the corresponding 
fractional contribution for the antiprotons. 



C. Systematic Uncertainties 

Systematic uncertainties on the hadron spectra are es- 
timated as in reference 0- Various sets of pt spectra 
and ratios of different particle types were made by vary- 
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TABLE II: Systematic errors in percent on particle yields in d+Au and p+p collisions. These values are independent of pr- 



d+Au p+p 

geometric acceptance correction 4 4 

track matching 9 8 

timing variations 5 5 

reconstruction efficiency correction 3 4 

energy loss correction 1 2 

trigger bias - 4 



g g ^ Fractional feed down correction f 

E protons 

"■^r antiprotons 

0.7 L L I Systematic error for protons 

0.6 - 




0.1^ 



"o 0.5 1 1.5 2 2.5 3 3.5 4 4.5 

Pj GeV/c 

FIG. 3: Fractional contributions of protons, f, (solid line) 
and antiprotons (dashed line), as a function of pr, from weak 
decays in all measured protons and antiprotons. Systematic 
error band (26%) for the protons is shown and discussed in the 
text. The same systematic error applies for the antiprotons. 



ing the cut parameters such as fiducial cut to check ac- 
ceptance corrections, track association (i.e. matching) 
windows, and PID cuts, from those used in the analysis. 
For each of these spectra and ratios the same changes in 
cuts were made in the Monte Carlo analysis. The un- 
certainties were evaluated by comparing fully corrected 
spectra and ratios from different cuts. The resulting un- 
certainties from each cut are given in Table II and added 
in quadrature to yield overall systematic uncertainties. 

Additional systematic uncertainties on the hadron 
spectra arise from time variations in the TOF timing 
(slat- by-slat and run- by-run variations), the small re- 
maining contamination by other species after matching 
and PID cuts, uncertainty in the corrections for track 
reconstruction efficiency, and uncertainty on particle en- 
ergy loss in the detector material. These uncertainties, 
which do not depend measurably on the hadron momen- 
tum, are listed in Table II. The sizable uncertainty in 
the matching cut is due to the non-Gaussian tails on the 
z-coordinate matching distributions. The track match- 
ing in this direction is limited by a relatively poor ver- 
tex resolution determined by the BBC in p+p and d+Au 
collisions due to the small multiplicities. The quoted 3%- 



4% uncertainty in the reconstruction efficiency correction 
represents the maximum local discrepancy between effi- 
ciencies measured with strict and loose track quality cuts. 

Uncertainties due to particle identification cuts are mo- 
mentum dependent. For protons and antiprotons, the 
identification uncertainty is 8% at low pr and decreases 
to 3% at high px- Kaons at low momentum have 10% 
PID uncertainty, decreasing to 3% at high px- For pions 
the uncertainty increases from 4% to 10 % with increas- 
ing pt- Kaon and proton uncertainties decrease with 
increasing px because energy loss and decay corrections 
become smaller. The pion uncertainties are dominated 
by the particle identification performance, which wors- 
ens with increasing px ■ 

The systematic error on the feed-down proton spec- 
trum is 26%, primarily due to uncertainty in the mea- 
sured A spectra and particle composition. The resulting 
systematic error on the final prompt proton and antipro- 
ton spectra is of the order of 10% in both p+p and d+Au. 
The systematic error on the proton to pion ratio is 12%, 
including the uncertainty on A/A. 

Systematic uncertainties on the d+Au nuclear modifi- 
cation factors mostly cancel as the p+p and d+Au data 
were collected immediately following one another, and 
detector performance was very similar. The overall sys- 
tematic error in the nuclear modification factor is due 
to uncertainties in the reconstruction efficiencies, fidu- 
cial volumes, and small run-by-run variations. It is ap- 
proximately 10%, independent of particle species and px- 
An additional d+Au scale uncertainty is shown as boxes 
around 1.0 in the figures; this is the quadrature sum 
of uncertainties on the p+p cross section of 9.6%, and 
the number of binary collisions in the each centrality bin 
(presented in Table I). 

The systematic error on the Au+Au nuclear modifi- 
cation factors is derived by propagating the systematic 
errors on p+p and Au+Au data 7] to the final ratio. The 
average systematic error for pions is approximately 15%, 
while for protons and antiprotons it is on the order of 
19%. The normalization uncertainty, as in d+Au, is the 
quadrature sum of uncertainties on the p+p cross section 
and the error on the number of binary collisions in the 
corresponding Au+Au centrality bin from reference 0- 
We note that for the most central Au+Au bin (0-5%), 
iVco;/=1065.4 and the uncertainty is ±105.3; in the most 
peripheral bin (60-92%), iVco;/=14.5 ± 4.0 . 



7 



>• 

73 10 


^ t. 

%. 


Z 1 

CM 


•A. 


10 




c 












■ 






lo- 





. • I • 



■ • 




PT(GeV/c) 



Pt (GeV/c) 



>. 

■D 

£10' 

10-^ 

£ 

'^^lo-^ 

c 

LU 

z lo" 

CM 

10 = 




Pt (GeV/c) 




FIG. 4: (color online) Invariant yields at mid-rapidity for 
positive and negative pions as a function of pr for various 
centrality classes in d+Au and p+p collisions. The error bars 
show statistical uncertainties only and are typically smaller 
than the data points. 



FIG. 6: (color online) Invariant yields at mid-rapidity for pro- 
tons and antiprotons as a function of pr for various centrality 
classes in d-|-Au and p-f p collisions. The error bars show sta- 
tistical uncertainties only and are typically smaller than the 
data points. 
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FIG. 5: (color online) Invariant yields at mid-rapidity for 
positive and negative kaons as a function of pr for various 
centrality classes in d-|-Au and p-l-p collisions. The error bars 
show statistical uncertainties only and are typically smaller 
than the data points. 



III. 



RESULTS 



A. Hadron spectra 

The fully corrected px distributions of tt, K, p, and p 
for the four d-f Au centrality bins and for p-l-p collisions 
are shown in Figures^ |S1 andlHl respectively. Pions show 
a power law spectral shape, while kaons and protons are 
exponential. 

In order to probe the hadron production mechanism, it 
is instructive to compare particle and anti-particle spec- 
tra. Figures [7| |H1 and |51 show the ratios of antiparticle 



to particle production as a function of pT in p-l-p, d-|-Au 
and, for comparison, central Au+Au collisions from refer- 
ence Q for TT, K and p, respectively. For all three hadron 
species the ratios are flat with pT- d-|-Au yield ratios 
are in good agreement with p-l-p collisions, and the ra- 
tios remain the same even in central Au-|-Au collisions. 
The production ratio of antiparticle to particle is 0.99 ± 
O.Ol(stat) ± 0.06(syst) for pions and 0.92 ± O.Ol(stat) 
± 0.07(syst) for kaons in both minimum bias d-f Au and 
p-|-p collisions. The antiproton to proton ratio is mea- 
sured to be 0.70 ± 0.01 (stat) ± 0.08(syst) in minimum 
bias d-l-Au collisions and 0.71 ± 0.01 (stat) ± 0.08(syst) 
in p-l-p collisions. All ratios are consistent within errors 
with values reported by PHOBOS and BRAHMS 

m. 



B. Nuclear Modification Factors 

The measurement of identified hadrons in both d-|-Au 
and p-|-p collisions allows study of the centrality depen- 
dence of the nuclear modification factor in d-|-Au. A stan- 
dard way to quantify nuclear medium effects on high pT 
particle production in nucleus-nucleus collisions is pro- 
vided by the nuclear modification factor. This is the ratio 
of the d-|-A invariant yields to the binary collision scaled 
p-|-p invariant yields: 



RdAiPr) 



{l/Nlf)<eNdAldydpT 
TdAud^all^JdydpT 



(5) 



where TdAu — {N cou) / crenel describes the nuclear geome- 
try, and d^a^^^il dydpT for p-l-p collisions is derived from 
the measured p+p cross section. {Ncoii) is the average 
number of inelastic nucleon-nucleon collisions determined 
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FIG. 7: (color online) Ratio of mid-rapidity spectra for n~ to 
tt"*" in d+Au, p+p and central Au+Au collisions. The error 
bars show statistical uncertainties only. 



FIG. 9: (color online) Ratio of mid-rapidity spectra for an- 
tiprotons to protons in d+Au, p+p and central Au+Au col- 
lisions 0). The error bars show statistical uncertainties only. 
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FIG. 8: (color online) Ratio of mid-rapidity spectra for K~ to 
in d+Au, p+p and central Au+Au collisions. The error 
bars show statistical uncertainties only. 



from simulation using the Glauber model as input, as 
described in section III Al N^"^ is the number of d+Au 
events in the relevant centrality class. 

Figure ^] shows RdA for pions, kaons and protons for 
minimum bias d+Au collisions. We observe a nuclear en- 
hancement in the production of hadrons with pt ^ 1-5 - 
2 GeV/c in d+Au collisions, compared to that in p+p. 
As was already suggested when comparing the enhance- 
ment for inclusive charged hadrons with that of neutral 
pions 33], there is a species dependence in the Cronin 
effect. The Cronin effect for charged pions is small, as 
was observed for neutral pions. The nuclear enhancement 
for protons and antiprotons is considerably larger. The 
kaon measurement has a more limited kinematic range, 
but the RdA is in agreement with that of the pions at 
comparable px- 

Figure ITTI shows RdA for pions, kaons and protons in 
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FIG. 10: (color online) Nuclear modification factor RdA for 
pions, kaons and protons in d+Au collisions for minimum 
bias events. The error bars represent the statistical errors. 
The box around 1.0 shows uncertainties in the p+p absolute 
cross section and in the calculation of NcoU- For the proton 
and antiproton RdA, the ~10% systematic uncertainty is also 
presented as boxes around the points. The systematic uncer- 
tainty on the pion and kaon RdA is similar but not shown in 
the picture for clarity. 



the four d+Au centrality bins. Peripheral d+Au colli- 
sions {{Ncoii) = 3.1 ± 0.3) do not show any modification 
of high momentum hadron production, compared to that 
in p+p collisions. At < 1 GeV/c, the nuclear mod- 
ification factor falls below 1.0. This is to be expected 
as soft particle production scales with the number of 
participating nucleons, not with the number of binary 
nucleon-nucleon collisions. More central collisions show 
increasing nuclear enhancement in both high pT pion and 
proton production. 

The bands in Fig. II II show a calculation of the Cronin 
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FIG. 11: (color online) Nuclear modification factor RdA for pions, kaons and protons in d+Au collisions in four centrality bins. 
The error bars represent the statistical errors. Boxes around 1.0 show uncertainties in the p+p absolute cross section and in the 
calculation of Ncoii- For the proton and antiproton RdA, the ~10% systematic uncertainty is also presented as boxes around 
the points. The systematic uncertainty on the pion and kaon RdA is similar but not shown in t he p icture for clarity. Solid 
bands show the calculation of the nuclear modification factors for pions by Accardi and Gyulassy [l^. 



effect for pions by Accardi and Gyulassy, using a pQCD 
model of multiple semi-hard collisions and taking geomet- 
rical shadowing into account jl^ . The agreement above 1 
GeV/c, where the calculation should be reliable, is very 
good for all four centrality bins. This agreement illus- 
trates that the multiple partonic scattering and nuclear 
shadowing alone can explain the observed Cronin effect 
and leaves very little room for gluon saturation effects in 
the nuclear initial state at mid-rapidity at RHIC [T^ . 



C. Centrality Dependence 

We further probe the effect of cold nuclear matter 
upon the hadron production using the number of col- 
lisions suffered by each projectile nucleon for the four 
centrality bins. Figure [T^ compares the centrality depen- 
dence of RdAu for pions and protons in two momentum 
bins. The modification factors are plotted as a function 
of = Ncou/Np^j.f, the number of collisions per partici- 
pating deuteron nucleon. The lower momentum bin, for 
0.6 < pt < 1.0 GeV/c, is chosen in the region where 
RdAu is less than 1.0, and hadron yields scale very nearly 
with the number of nucleons participating in the colli- 
sion, rather than with the number of binary collisions. 
As expected, RdAu decreases with v in this pr range, 
with negligible difference between pions and protons. In 
the higher px bin, RdAu increases with the number of col- 
lisions, with a notably larger rate of increase for baryons 



than for mesons. Though the RdAu values for higher pT 
hadrons appear to flatten with increasing centrality, the 
uncertainties are too large to allow a definitive conclusion 
about saturation with the number of collisions suffered 
by each participant nucleon . 
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FIG. 12: (color online) Integrated RdA for pions and pro- 
tons in two momentum bins as a function of the number of 
collisions suffered by the deuteron participant u. Error bars 
indicate the quadrature sum of statistical errors and uncer- 
tainties on the number of collisions bin-by-bin. The solid box 
on the left shows the magnitude of the centrality independent 
uncertainties. 
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FIG. 13: (color online) Ratio of the proton over the pion 
nuclear modification in d+Au collisions, for central and pe- 
ripheral events. Error bars indicate statistical errors only. 



D. Cronin Effect For Baryons 

Figure [TSI shows the ratio of the nuclear modification 
factors observed for protons and antiprotons to that for 
pions in the most central and the most peripheral d+Au 
collisions. The enhancement for protons is stronger, by 
approximately 30%-50% in the most central collisions; 
this was also reported by STAR ^] . The increasing dif- 
ference between baryons and mesons for more central col- 
lisions indicates that the baryon production mechanism 
appears to depend upon the surrounding nuclear medium 
already in d-|-Au collisions. We note, however, that the 
species dependence of the Cronin effect in d-|-Au colli- 
sions is much smaller than the factor of «3 enhancement 
of protons in central Au-|-Au collisions, as can be seen 
by comparing the nuclear modification factors for pions 
and protons in central Au+Au collisions shown below in 
section FlVAl 

By converting the nuclear modification factors to per- 
beam-nucleon cross section ratios and vice versa, it is 
possible to compare to measurements of the Cronin ef- 
fect at other energies. Figure El shows nuclear modifi- 
cation factors from this work compared to those derived 
from a-factors measured at lower energies 0] , in a similar 
manner to that used in |^ to calculate per-nucleon cross 
section ratios. In order to make the transition to nuclear 
modification factors from per-nucleon cross section ratios 
we assume that cr^+A = 2 x ap-^-A and Cp+d = 2 x (Tp+p 
at the low energy of interest, which is a very reasonable 
approximation for all particle species j^j- 

The observed species dependence of the enhancement 
is similar to that measured in lower energy collisions Q. 
The magnitude of the enhancement for pions at > 3 
GeV/c is larger at ^ = 27.4 GeV than at 200 GeV. Pro- 
tons and antiprotons are also more enhanced at the lower 
beam energy: a factor of 3.5 at ft! 4 GeV, as compared 
with a factor of 2 at ^/s = 200 GeV. This energy depen- 
dence of the Cronin effect for pions has been interpreted 
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FIG. 14: Nuclear modification factors for charged pions, neu- 
tral pions '3^, and protons and antiprotons from minimum 
bias d-fAu from this work (shown by open symbols) compared 
to nuclear modification factors calculated from the per-beam 
nucleon cross sections reported for ^/s — 27 A GeV p -|- A 
collisions 2] (closed symbols and continuous bands) 



as evidence for a different production mechanism for high 
Pt hadrons at RHIC compared to lower energies j^ll . In 
this model, high px hadrons are produced incoherently 
on different nucleons at low energy, while in higher en- 
ergy collisions the production amplitudes can interfere 
because the process of gluon radiation is long compared 
to the binary collision time. Coherent radiation from dif- 
ferent nucleons is subject to Landau-Pomeranchuk sup- 
pression. However, the difference between baryon and 
meson Cronin effect is not predicted by this model. 

IV. DISCUSSION 

Traditional explanations of the Cronin effect all involve 
multiple scattering of incoming partons that lead to an 
enhancement at intermediate pt 12]. There are vari- 
ous theoretical models of the multiple scattering, which 
predict somewhat different dependence upon the number 
of scattering centers. The observed centrality or v de- 
pendence for pions is well-reproduced by semi-hard ini- 
tial state scattering as shown in Fig. ^2 see also 
0, 0, 0|. The models include initial state multiple 
scattering as well as geometrical shadowing. However, 
none of these models would predict a species dependent 
Cronin effect, as initial state parton scattering precedes 
fragmentation into the different hadronic species. The 
markedly larger Cronin effect for protons and antipro- 
tons requires processes in addition to initial state multi- 
ple scattering in baryon production at moderate trans- 
verse momenta. 

Recently, Hwa and collaborators 0| have shown an 
alternative explanation of the Cronin effect, attributed 
to the recombination of shower quarks with those from 
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FIG. 15: (color online) The ratio of feed-down corrected pro- 
tons to TT^ and antiprotons to 7r~ in minimum bias p-l-p and 
d-f An compared to peripheral and central Au+Au collisions. 
Statistical error bars are shown. 



the medium in d+Au collisions. Such models do predict 
a larger Cronin effect for protons than pions, and may be 
justified by the rather short formation times of = 2-4 
GeV/c protons. According to the uncertainty principle, 
the formation time in the rest frame of the hadron can be 
related to the hadron size Rh- In the laboratory frame, 
the formation time of a hadron with mass rn/i and energy 
Eh is given by r/ « i?/^^. For & pr = 2.5 GeV/c pion, 
the formation time is 9-18 fm/c (for Rh — 0.5 - 1.0 fm), 
well outside the collision region. However, for pT = 2.5 
GeV/c protons, the corresponding formation time is only 
2.7 fm/c, suggesting that the hadronization process may 
well begin in or near the nuclear medium. 



A. Comparison to Au+Au Collisions 
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FIG. 16: (color online) Nuclear modification factors for pi- 
ons, comparing central and peripheral Au+Au collisions to 
central d+Au. It should be noted that the number of binary 
nucleon-nucleon collisions in peripheral Au+Au and central 
d+Au is very similar. Solid bars on the left indicate normal- 
ization uncertainties in the p+p absolute cross section and in 
the calculation of Ncoii for the three systems. Error bars indi- 
cate statistical errors only, while for the most central Au+Au 
case the systematic errors, discussed in the text, are shown as 
boxes around the points. 
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The proton to pion ratio from minimum bias p+p and 
minimum bias d+Au are compared to each other and 
to central and peripheral Au+Au collisions in Fig. ^1 
As noted above, protons and antiprotons are feed-down 
corrected in each system. 

The p/tt ratio in d+Au is very similar to that in periph- 
eral Au+Au collisions, and lies slightly above the p+p 
ratio. The p/ir ratio in central Au+Au collisions is, how- 
ever, much larger. The difference between the ratio in 
d+Au and central Au+Au clearly indicates that baryon 
yield enhancement is not simply an effect of sampling a 
large nucleus in the initial state. The large enhancement 
requires the presence of a substantial volume of nuclear 
medium. 

Figures lTCl 1 1 7l andl 1 8lcompare the nuclear modification 
factors for pions, kaons and (anti)protons in Au+Au and 
d+Au collisions. The p+p data from this work allow. 



FIG. 17: (color online) Nuclear modification factors for kaons, 
comparing central and peripheral Au+Au collisions to central 
d+Au. Statistical and systematic errors are represented as in 

Fig.rrsi 



for the first time, the calculation of nuclear modification 
factors in Au+Au by PHENIX. The Au+Au data are 
taken from reference J\ . It should be noted that common 
fluctuations between RdAu and RauAu in Figures [T7I and 
1181 arise because the p+p denominators are common. 

Central and peripheral Au+Au collisions are compared 
to central d+Au collisions, which have a similar number 
of binary collisions as the peripheral Au+Au sample. Pi- 
ons show a much lower Raa at high p^ in central than in 
peripheral Au+Au collisions, as expected from the large 
energy loss suffered by the partons in central collisions. 
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tion factor in d+Au for protons shows a larger Cronin 
effect than that for pions, and the difference increases 
with coUision centrahty. This difference was seen, but 
never fuUy understood, in lower energy collisions, and is 
not large enough to account for the abundance of protons 
in central Au+Au collisions. The difference between pi- 
ons and protons does, however, indicate that the Cronin 
effect is not simply due to multiple scattering of the in- 
coming partons. Raa for protons and antiprotons con- 
firms previous observations that the production of high 
Pt baryons in Au-|-Au scales with the number of binary 
nucleon-nucleon collisions, but the baryon yield per col- 
lision in Au+Au exceeds that in p+p. 



FIG. 18: (color online) Nuclear modification factors for 
protons and antiprotons, comparing central and peripheral 
Au+Au collisions to central d+Au. Statistical and system- 
atic errors are represented as in Fig. ^| 

The nuclear modification factor rises faster with px in 
d+Au than in peripheral Au+Au, despite the compara- 
ble number of binary collisions. As Au+Au involves a 
second Au nucleus, shadowing effects can be expected to 
be larger, reducing the observed Cronin effect. 

The proton and antiproton nuclear modification fac- 
tors show a quite different trend. The Cronin effect, 
larger than 1.0 at higher pT values, is independent of 
centrality in Au+Au collisions. This feature was already 
observed as binary collision scaling of proton and antipro- 
ton production in the central/peripheral collision yield 
ratios j^. The Cronin effect in d+Au is at least as large 
as in peripheral Au+Au. The difference indicates that 
baryon production must involve a complex interplay of 
processes in addition to initial state nucleon-nucleon col- 
lisions. 



V. CONCLUSION 

We have presented the centrality and species depen- 
dence of identified particle spectra in d+Au collisions, 
including the dependence of the nuclear modification fac- 
tor upon the number of collisions per participant nucleon. 
We also presented the first measurement of Raa for pi- 
ons, kaons, protons and antiprotons in Au+Au collisions. 

The Cronin effect for charged pions is small, but non- 
zero. The proton to pion ratio in d+Au is similar to 
that in peripheral Au+Au, while the corresponding ra- 
tio in p+p is somewhat lower. The nuclear modifica- 
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